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The present report covers essentially the pre-launch phase of the 
TELLUS Project. 

In this period, algorithms for the interpretation of remotely sensed 
surfaa© temperatures were developed (€.2.1.3). These are the TERGHA 
model (6.2. 1.3.2) for the calculation of evapo transpiration and soil 
moisture of vegetated areas and the TELL-US model (6. 2. 1.3. 3) for the 
evaluation of evaporation and soil moisture of bare soil or scarcely 
vegetated areas. Additional methods for the calculation of regional 
evapotranspiration were also evaluated (6.2. 1.3.4) . 

Experience in the interpretation of thermal data was obtained in 
flight campaigns in the Netherl^ds, France and the United Kingdom 
and in long term field measurements in France and Italy (6.2.2). 

These caunpaigns permitted the test and application of the interpreta- 
tion algorithms and permitted the development of data processing 
routines . 

The use of HCMM data for studies pi a regional heat budget was 
simulated from aircraft data, Likewise, the identification of frost 
prone areas and the influence of wind on remotely sensed surface 
temperatures have been investigated (6.3.1). 

The HCMM data products received were generally of good quality (5.1.3). 
Only a. visual interpretation regarding geomorpho logic and topographic 
features has been carried out so far (6. 2. 1.2). 


1 . INTRODUCTION 


The aain objective of the TELLUS Project is the development and pre- 
sentation of remotely sensed data on earth resources. This will permit 
the Commission of the European Communities (EC) to comparatively evalu- 
ate the potential and benefits of air-and-space observation systems 
for obtaining information required for viurious sectors of its common 
policy. 

An attempt will be made to learn whether HCMM data/ together with 
other existing aircraft and satellite data had a significant advantage 
over conventional methods for providing the Commission with this in- 
formation. In this sense, the TELLUS Project should give a practical 
demonstration on a European basis. Thus it is being developed as a pilot 
experiment aimed to cover a range of significant European conditions 
regarding agriculture, hydrology, geology and the environment. 

Furthermore, set up as a collaboration of several national Institutes 
and Organizations working together under the leadership of the EC's 
Joint Research Centre (JRC) Ispra Establishment, the TELLUS Project 
permits the concentration of diversified competences and support faci- 
lities for the comparative study of various typical European tet't sites. 
This form of collaboration was chosen both for reasons of efficiency and 
for promoting Remote Sensing methods within the European framework. 

The objectives of the proposed programme , while partially fulfilling a 
number of specific EC needs in the policy of agriculture and the environ- 
ment, also correspond to ,the specific research objectives of the parti- 
cipating Institutes. 


2. FRAMEWORK OF COINVESTIGATORS 


The list of the Organizations and Institutes participating in the TELLUS 
Investigation (HCM-025) is reported in Table 1, together with the names 
of Co- Investigators. 

The Co- Investigators are divided into six groups: 

SC General Directorates and JRC (Ispra Establisiunent) 

Benelux Institutes 
British Institutes 
French Institutes 
German Institutes 
Italian Institutes. 
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TABLE 1 : .LIST OF ORGANISATIONS AND INSTITUTES PARTICIPATING IN THE 
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TELLUS PROJECT (HCM^25) 
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TEST-SITE 

COORDINATORS 

(TSC'SI 

CO-INVESTIGATORS 

(COI'S) 

ORGANIZATIONS AND INSTITUTES Q^IQINAL •’ 




Oiraciarataaoltlia Europian Communitlaa 





; Oiraeioiiii Qmaral for Afrleullura, Bmaali 

OG-VI 




• OIrteionta Ganaral for Rittarcb, Scianea and Eduation, Bnjaaaii 

D5-XII 




• Joint Rawareb Contra, iRiri 

JRC 




• EEC Oalafation, Widiinpan 

OG-f/WO 




National Orftaimiant and Inninilaf 

1 

* 

FtddM A.A. 


Stlfium 

1 

1 



O'Hoori J.M. 

• Fakultaitdaf landbouvm«ttaniehaopan, Kath. Univaraittit, Lauvan 

KUL/LBB i 




Oinmirk 

1 



Ertirdt F.£,‘ 

^ ■ Inititutft for flkologisk Botanik, Univanjtv of Cooanhagan 

UK/IQB 1 

i 


Solllot Ch. 


fnna 




Geillot Cli. 

• Sarvica da TilMdtanlon, INRA Vwaillat 

INRA/ST 



Birloy E. 

■ Station d'Amalioration dav Flantai, INRA Rannat ' 

INRA/SAP 



Etfuin B. 

• Station da Bioelimatologia, INRA Monrfivat 

INRA/CRASE/SE 



Firritr A. 

• Station da Bioeiimatolagfa, INRA Varaillai 

INRA/SB 



Sickir F. 

• Groupa da Raeharehaa an Tiladataetion Radiomdtrioua, 
Univartita da Stnubourg 

ULP/GRTfi 


Sohnil HJ. 


Stntmttf 


^ , 


Gonmmn H., Wtiiehtt W. 

• Saognpnltchta Inaiitut dar Univartitat, Fratburg 

UF/GI 



BBhnii H-i. 

. • Innitut fur Phviikailteha Waltraumforichung, Freiburg 

ipyy 


R. van darPloag 

■ Inatitut fur Bai.ankunda und Waldrnahrung 

“UG/IBW. _ 


Ctliindro A. 


nvy 



Boiriallo L. 

• Camro Snidi Aoplieazloni T acnologia Avanuta, Sari 

eSATA 



Cavtm U 

• litttuto di Agronomia Ganaraie a Coltivazioni Erbaeaa, 
Univariiti di Bologna 

UBO/IAGCE 



Miltlla A. 

• Istituto di Agronomia a Coltivazioni Arboraa, Univarsita di Saaari 

USS/IACA 



Caliindro A.. Faeucci G, 

■ Istituto jfi Agronomia a Coitivazionl Erbaeaa, Univarsita di Barf 

USA/IACE 

REPRODUCIBILITY OF 

' Istituto di Agronomia Ganaraia a Coltivazioni Erbaeaa, 
Univarsiti di Firanza 

UFI/IAGCE ' 

ORIGINAL PAGE IS POORp,«f. 

• Istituto di Fisica, Univrsita di Bari 

U BA/IF 



Martalonyo 3. 

• Istituto di Gaologia Aoolicaia, CNR Padova 

IGA i 



Fiarntiiiprina A. 

• Istituto di Minaralogia a Gaologia, Univarsita di Sassari 

USS/IMG ! 



Ueni Q.M. 

• Istituto par la Gaofitica dalla Litotfara, CNR Milano 

IGL i 



Tombas L 

• Inituto Sparimantala par la Nutrizlona della Pianta, Roma 

ISNP 1 



Roant E. 

■ Uffieio Cantrale di Ecologia Agrarla, Roma 

UCEA ‘ 

* 

Fiddti R.A. 


Ntthirknds 


y 


Faddai R.A. 

Innituut voor Cuttuurtachniak an Wnarbuidiouding, Waganingan 

ICW 

■ 


WluaJ. 

• Koninklijk Nadariands Metaorolagiacb Instituut, Oa Silt 

KNMI 


Me Culloch J.G. 


Unittd Kingdom 




Kirkby MJ. 

■ Department of Guograony, University of Leads 

UL/OG 



Savigaar R.A.G. 

■ Dapartmant of Gsograohy, University of Reading 

Ufl/DG 



Me Cdlloeh J.G. 

■ Institute of Hydrology, Wallingford 

IHW 
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3 . TEST SITES 


The TBLLUS test sites are reported in Table 2. Geographic as well as 
technical considerations led to regroupe them in 5 groups according 
to the Institutes involved and to designate for each group one Co- 
Investigator as test site Coordinator. Names of test-site Coordinators 
are reported in Table 1. 

Pig. 1 shows the location of the TELLUS test sites within the European 
coverage area of the HCMM receiving stations of MADRID (NASA) and 
LANNION (ESA). The "diamond" on the right hand side of Fig.l encloses 
twelve hour night/day coverage pattern due to the intersection of the 
corresponding satellite or-;ites. 


4. OBJECTIVES OF THE INVESTIGATION 

According to the interest of the EC's Commission and the proposals 
made by the national Institutes, three main thematic lines of investi- 
gation have been singled out: 

1. Evaluation of evapotrauispiration and moisture content of bare soils 
and soils covered by vegetation 

2. Study of the interactions between natural phenomena and mesoscale 
heat budget 

3. Man-made changes and their impact on regional heat budgets. 

Tables 3 - S summarize the research subjects and objectives of the 
TELLUS Co-Investigators . 

Estimates of soil moisture content and of evapotranspiration can lead 
to considerable cost savings, by a more effective planning of water 
supply, irrigation and drainage systems. 

The knowledge of soil moisture content is important for both dry areas 
and for moisture satiirated zones. In the first case, lack of water in 
soils narrows the range of profitable crops, causes crops stress and 
reduces crop yields. 

In the second case, saturation of soil is important both from an agri- 
cultural point of view and for monitoring flooding. Accelerated erosion 
is also influenced by the distribution of areas of high soil moisture 
and the depth of the ground water table. 

The heat budget is a relevant environmental characteristic due to its 
impact on vegetation and man. Agricultural productivity may be strongly 
affected by variations in the regional heat budget due to both man-made 
and natural (i.e. meteorologic) causes. 

Urban and industrial areas (factories, power stations) can significantly 
affect the heat balance of water bodies and surrounding farmland and 
woodland. The evaluation of the heat budget on a regional scale with 
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CAPACITY MAPPING MISSION (HCMM) 12-HOUR NIGHT/ DAY COVERAGE 
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Table 4 Man-made Changes and their Impact on Regional Heat Budget 





on Regional Heat Budget 
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mixed rural and industrial activities is one of tho important steps in 
the conservation of environmental quality. 

Evaluating and monitoring the tendency to frost of agricultural areas 
and correlating between thermal, hydrological and topographical cha- 
racteristics of band surfaces and the occurrence of hail stoinns are two 
among the subjects for investigation where not man, but nature influences 
the regional heat budget. 

3. DATA SECSyriON hNO DISTRIBUTION 

HCMM data reaching JRC can be divided into two groups: 

1. Data as stipulated in the Provisions for Participation i.e. test or 
preliminary data, standing order data (SO) and retrospective order 
data (RO) 

2. Quick- look (QL) data as s*tated in the Freiburg Meeting of WG 3 
(Dec. 1973). 

5 . 1 Data according to Provisions for Participation 

JRS has received only test data from MASA-GSPC until the end of 1973 in the 
form of positive prints and transparencies 24x24 cm, DIE/DVIS and NIR and 
a number of corresponding tapes. 

Delivery of standing order format data in the form of negative 24x24 
transparencies started at the beginning of 1979. A complete list of 
SO data received is given in Table 6 . 

No RO data has been received up to now. 

5.1.1 Data transmission 

A flow-chart of the data transmission is presented in Pig. 2 . kCMM data 
are sent by regular mail from GFSC to the SC Washington Delegation. Prom 
there they are forwarded by diplomatic pouch via Brussels to JRC Ispra. 
Travel time frcm GSPC to JRC did not exceed 16 days. 

Dispatching of the data to the various test site coordinators is carried 
out at Ispra as follows: 

In order to localize the test sites, SO transparencies are superposed on 
a transparent map of Europe carrying the geographic coordinates and 
the contours of the test sites. 

Problems have been arisen to crrectly localize images with heavy cloud 
cover when erroneous coordinates were indicated on the images. 


TABLE 6 ; 


PHority Proctsting List (PPL) and ihipmanti rtc( jvad at JRC 


AL 

OaM 

Tatt-iitas 

Track 

AL 

Data 
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14 
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B : Benelux 

0 ; Germany 
F : France 
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FIG. 2 : HCMM DATA DISTRIBUTION FROM GSFC 
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5.1.2 Data received 

The first SO received was processed on January 25, 1979, t.he last one 
on June, 1979. 

The data recived correspond to 49 orbits relevant to the investigation and 
are distributed as listed in Table 6 . in this table, the "presence" of a 
test site group indicates that at least a part of one test site of the 
group was covered by the track, regardless of cloudiness. 

5.1.3 Data quality 

reproducibility of the 

Data quality is ummed up in Table 7 , ORIGINAL PAGF IS POOR 



test images 

SO images 

Last SO images 

Scanning 

good 

good 

good, 

intermittent poor 

Optical 

quality 

good but 
low contrast 

good 

irregular 

Dens i tome trie 
quality 

low 

very good 

good when 
scanning is good 

Cloudcover 

irregular 

intermittent 

very irregular 


Table 7 : Data quality 


5.1.4 General remarks 

The end of January 1979 GSFC asked JRC to produce a priority processing 
list (PPL) because of the. long delay in SO delivery. 

This list (Table 3 ) was established according to the requests of the 
co-investigators and was sent to GSFC in February 1979. 

Up to now only 20 out of 80 priority requests have been received. Inquiries 
as to the criteria used in processing images received no reply. 

The PPL was also communicated to Lannion with a request for retrospective 
quick-looks. This request could not be fulfilled as the station is at 
present fully occupied with the production of 1979 QL. 


TABLE 3 ; 


Prtianct of Europaan groups of tast litai on tha S.O. I'magas transmittad to JRC Ispra, than dispatehad 
to tha various co*invcstigators 
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+ 
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+ 
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07.02 

+ 
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+ 

+ 

+ 
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+ 
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+ 
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+ 


+ 

+ 
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+ 





5 
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+ 

+ 
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+ 

+ 

+ 

+ 
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07.05 

+ 

+ 

+ 

+ 
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+ 

+ 

+ 
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+ 

+ 
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6 
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+ 

■ + 
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+ 
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12.29 
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3 
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+ 

+ 

11 
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+ 

+ 

+ 
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+ 


+ 

+ 
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+ 
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+ 

14 

05.11 


+ 



+ 






08.15 
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08.31 











10.18 
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1 
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12.15 








+ 

X 


12.21 







+ 



IS 

07.31 




I 
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09.01 






+ 


4- 

+ 


11.04 








+ 

+ 
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+ 


t 

X 

T 


12.22 








+ 



Consatiutiva 

night/day 


10 IS 7 12 




TOTAL 


19 13 13 13 4 


1 


5.2 Qulck-^ looks 


After the signature of an agreement beta^een SSA and the French government 
in May 1979, CMS Lannion staurted to produce negative transparent QL. They 
are dispatched regularly at fortnightly intervals from CMS Lannion di- 
rectly to the various test site coordinators. A flow diagram of the 
distribution is shown in Fig. 3 . 

The quality of the data is good. QL are produced irrespective of cloud 
conditions . 


6. ACTIVITY PERFORMED A^^D RESULTS OBTAINED 


SS-S.*" 


6.1 Generalities 

In 1976 the TELLUS Project was decided to be conducted in two phases: 

1) pre-launch, 2) post- launch and continuing investigations. 

However, owing to the long delay in receiving HCMM data (see 5.) and 
the time-shift of .the official start of the HCM-02S Investigation, the 
present report covers the "TELLUS First Development Phase" which goes 
somewhat beyond the launch date. 

During this 'phase research was aimed at providing the different co-in- 
vestigators with the necessary tools for solving the problems connected 
with the future HCMM data interpretation. 

Depending on the principal subjects of investigation a number of prob- 
lems were considered, depending on -their importance mainly in respect 
to the methodological aspects and on the possibility of tackling -them 
with the available technical means. 

The principal effort was devoted to determine the me-thodologies in 
ground measurements, airborne survey auid in data processing and inter- 
pretation. Priority was given to setting up ma-thematioal models and 
visual interpretation techniques of satellite data and to problems of 
data handling such as geometric corrections and superposition of night- 
and-day aircraft data. 

6.2. Evaporation and Soil Moisture 

6.2.1. The approach adopted 

6 . 2 . 1 . 1 . Me-thodological aspects 

A two-way approach has been adopted; 

1. A first approach consists of the visual interpretation of the HCMM- 
images (originals or transformed) . 

The method may be described as an analog-deductive interpretation pro- 
cedure, final products being hydro-geological maps. 










2. The second one involves a mathematical modelling procedure. 

It employs digital models the parameters of which are obtained from 
air-and-space platforms (aircraft and/or HCMM satellite) and ground- 
based instrumentation. In fact, one method of estimating evapotrans- 
piration and soil moisture content is to use a model of the energy 
balance of the soil surface and of the atmospheric boundary layer. 
Such a model may be written in a form that most of the input data 
used are routinely observed and it may be adapted so that remotely- 
sensed measurements may be included. 

Thus it may be considered a mathematical-experimental interpretation 
approach. Final products are thematic maps (evapotranspiration, soil 
moisture, etc.). 

This approach was chosen as the most suitable one for our investi- 
gation. In fact, synoptic estimation by remote sensing of evapo- 
transpiration and soil moisture is an important research subject 
of quite general interest not only among TELLUS Co-Investigators 
but also more extensively in Europe. 

Finding empirical correlation only between surface temperature and soil 
moisture and ET would perhaps suffice on a. local basis but- would not 
permit a comparative study of results obtained under the various 
European conditions. 

Therefore it was decided to gain a deeper understanding of the processes 
relating evapotranspiration and soil moisture with other ground para- 
meters and remotrsly sensed signals. The approach chosen was: 

1 . to obtain an exhaustive phenomenologic description by a mathematical 
model , 

2. to test the model for a number of field conditions. 

In this way, after confirmation that the mathematical expressions cor- 
rectly describe the processes involved, one can verify the importance 
of the various parameters and try to simplify the model accordingly. 

One may also attempt to extend the model from the punctual (local) 
condition to a more extended (regional) one. A full sensitivity analysis 
may also be applied to the model in order to identify the parameters 
which may be neglected, in relation to the reciprocal importance, the 
accuracy requested by the output and the nominal error of the remotely 
sensed signals (aircraft and/or satellite) . 

At this point, the different field tested mathematical expressions may 
be compared with the empirical relations obtained for the same field 
conditions employing regression analysis or similar techniques. This 
will be in favour of a better understanding of the processes involved. 
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We are conscious that the agricultural land structure in Europe is rather 
fragmented and small as compared to th^ HCMM resolution of 0.6x0. 6 km. 

The complex pattern of geomorphology and land use on the Eiuropean 
continent poses the problem of the meaning of HCMM pixels under these 
conditions. The thermal signal from a surface of 0.6x0. 6 km will in many 
cases be the result of the thermal emission from a number of small patches, 
each of them having different emission properties. 

An attempt is being made of "scaling up" the partial informations of the 
small pixels composing a HCMM pixel to the overall 0.6x0. 6 km information. 
This technique: is based on the comparison of groups of pixels of an air- 
craft flight taken at a given altitude with groups of pixels of an air- 
craft flight taken at a higher altitude over the same area. Through 
successive comparisons of increasing area between pixels corresponding 
to flights at higher altitudes one can arrive to' the final comparison 
with a satellite pixel. 

6. 2. 1.2 Visual interprets tion 

Visual interpretation of HCMM images was carried out using current 
methods of photointarpretation. A flowchart of the method is presented 
in Pig. 4 . 

The instrumentation employed consisted of stereo plotter, densitometer, 
analogs image analyser and an additive colour viewer. Maps, LANDSAT 
and Skylab images were used in the process and their relative effec- 
tiveness in identifying geomorphologic and topographic features was 
compared. 

A preliminary analysis of the first thermal test image from Cheasapeak 
Bay (May 11, 1978, 14.03 h) showed the superiority of the HCMM thermal 
image over LANDSAT images in detecting alluvial plains, marshes, low 
lying areas and relict river beds. LANDSAT was superior in identifying 
main water courses and water pollution plumes. A TSLLUS Newsletter on 
this work is in preparation. 

A further physiographic analysis was carried out for the Po Valley 
using VIS, DIR and NIR images as well as images from LANDSAT (bands 4, 

5 and 7) and Skylab (band 7) (Barisano, 1) . Regarding the hydro- 
logic aspects, only the DIR image permitted to identify the upwelling 
zones of aquifers (moisture in depth) , while the NIR image allowed 
localization of surface moisture. To a lesser decree the NIR image 
could be used to differentiate rock formations . 

A comparison of the usefulness of HC!^M, LANDSAT and Skylab data for 
identifying geomorphologic and topographic features in the Po Valley 
is given in Table. 9 . 
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analogic or 
digital 
processing 


interpretative keys 


COMPARISON 


LANOSAT images original 
or suitably processed 
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Figure 4 Flowchart of the visual interpretation process 


One should note ..the superiority of KCMM thermal images in identifying 
ground water, alluvial plains, moraines and ancient river beds. No 
general value should be attributed to the classification in Table 9. 
as it is limited to the case actually treated. 


TaUsle 9 .Comparative evaluation of HCNM, IiANDSAT and Skylab images 
for identifying geomorphologio and topographic features in 
the Po Valley 
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S . 2 . 1 . 3 Mathemitical modelling 
6 . 2 . 1 . 3 . 1 General description 

Two models have been developed up to now. The first one, the "TERGRA,” 
modal, developed by Soer (Soer, 2) is for use in grass land- type areas and, 
using one measurement of the soil surface temperature, gives estimates 
of -Jie daily evaporation. The other, the "TELL-US" model, developed 
by EARS (EARS, 3) , is for use on bare soil, and estimates both daily 
evaporation and thermal inertia using measurements of both day-time 
and night-time surface temperatures. 

Both models 2 u:e layered models, which require known boundary conditions 
at a given height in the atmospheric boundary layer and at a given 
depth in the soil in order to solve the equation of the energy balance 
of the crop-or-soil surface. The models are based on the equation of 
the heat diffusion in the soil 

^ i ^ REPRODUCIBILITY OF THE 

ORIGINAL PAGE IS POOR 

The atmosphere and the soil regime are coupled by the continuity re- 
quirement: 

R +G+K+LE-0 /20-1/ 

n 

which is the energy balance of the soil or crop surface. 

In /7-1/ R is the net radiation flux, G the soil heat flux, H the 
sensible heat flux in the atmospheric boundairy layer and LE the latent 
heat flux. The equation may be solved numerically for each level z ' to 
give a temperature profile using a finite element or finite difference 
method, tadting the soil or crop surface temperature as the upper 
boundary condition, and some knowp temperature at a given depth in the 
soil as a lower boundary condition. 

The net radiation flux is the difference between the incoming and outgoing 
short- and long-wave radiation fluxes 

R a (1 - a ) R -i- (1 - a-) R- - £ (7T^ /20-2/ 

n 3 s -C 4 o c 

where R is the incoming short wave and R^the incoming long wave radia- 
tion, a and <2^ the short and long-wave soil or crop reflection coeffi- 
cients s the soil or croo emission coefficient, (7 the Stefan 
Boltzmam constant, and T^. the soil-or-crop surface temperature. For long 
wave radiation- £ = i - 0.-f , so ecuation /20-2/ may be written 

Q i¥ 

4 

R a(l-a) R + srR, - (7T) /20-3/ 

n s s c C c 






The soil heat flux G through any z > z' is related to the temperature 
gradient in the soil (dT/dz) and the thermal cunductivity A , both 
evaluated at z * : 


G - 



/21-4/ 


The sensible heat flux may be estimated using a transport equation de- 
pending on the air temperatxire T^, the soil or crop surface, T^, the 
aerodynamic resistance, r^ and the air density Q and the specific 
heat Cp: 


H 



(T - T ) 
a c 


r 

a 


/21-5/ 


The latent heat flux may also be expressed as a transport equation: 

5c (e - e ) 
a s 


LE 


r 


r + r 
a s 


/ 21 - 6 / 


where y is the psychometric constant, the water vapour pressmres 

in the atmosphere layer, e the saturated water vapour pressure at 

temperature T , and r the stomatal diffusion resistance to water 
c s 

vapour transport. The quantities e and e may be estimated from T 
and Tjj. Quantity will be discussed in the following section on the 
TERGRA model. 

By comparison of the above terms it may be seen that the components of 
the energy balamce may be expressed as functions of the air temperature, 
the temperature at a given depth in the soil and the soil surface 
temperature, as well as of the incoming radiation and the aerodynamic 
and stomatal resistances. As said in 6. 2. 1.1, most of these variables 
are feeing measured or may be estimated from routine observations. The 
exception is the crop or soil surface temperature, T^, which is not 
routinely measured except at a few sites. 

This measurement has to be supplied by KCMM, by an airborne scanner 
or by ground based infrared thermometers. 


- Application of the model 

The initial estimate of the surface temperature T^ will probably not 
satisfy the continuity requirement of the energy balance equation, be- 
cause of the simplifications necessary to estimate a number of para- 
meters such as r^ or r^. Therefore, it is necessary to use an iterative 
numerical optimization technique to ensure that the condition of the 
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energy balance equation is satisfied. This procedure may be adopted to 
give estimates of the surface temperature T corresponding to every 
time step of a, set of routine measurements of air temperature, wet 
bulb temperature, wind speed and incoming short and long wave radiation. 

As said before, the soil or crop surface temperature measurements are 
not made routinely. As thermal emission measurements they are made 
instant 2 tneously by the radiometer mounted on the aircraft or on the 
HCMM spacecraft. Zt is possible to compare them with those estimated 
from the model and then to adjust the parameters of the model until 
the estimated surface temperatures are the seune as the measured surface 
temperatures at the same time. These parameters such as the aerodynamic 
resistance r^^ and the thermal heat capacity of the soil c^ may then be 
used to give more accurate estimates of daily evapotranspiration or soil 
moisture. 


6.2. 1.3.2. The TERGSA model 
- Specific motivation 
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In Western Europe, potential evapotremspiration of crops generally exceeds 
effective rainfall in spring and early summer. Under these cireuastemces , 
when the available moisture in the root zone is exhausted and capillary 
rise from the groundwater table does not succeed in supplying sufficient 
water to achieve the potential evapotranspiration rate, evapotranspira- 
tion will be reduced. This results in a lower crop growth and crop pro- 
duction. In such conditions, simultaneously with the decrease in evapo- 
transpiration, the temperature of the crop surface will increase. As 
distinct physical relation between the surface temperature and the 
evapotranspiration rate of a crop does exist, it is possible to calculate 
actual evapotramspiration from the crop surface temperature. Subsequent- 
ly, the moisture condition of the soil” cam be estimated from the evapo- 
transpiration rate. 

Up to now, owing to the high cost of measuring actual evapotranspi.- 
ration, conventional evapotranspiration measurements were limited to 
small hydrologically uniform areas. The use of remotely sensed crop 
surface temperatures eliminates the need of intensive field measure- 
ments and offers the possibility to obtain regional evapotranspiration 
estimates, making water balance analysis for large areas more feasible. 
For specific regions, crop production as correlated with the actual crop 
transpiration rate could be estimated. The extremely dry summer of 1976 
affecting large parts of Western Europe emphasized once more the impor- 
tance of having such large-area information. 

- Model description 


The TSEGHA model was developed as an aid for the interpretation of infra- 
red line scanning images of cropped surfaces, with particular emphasis on 


Matar and heat transport in the atmosphere are mainly passive transport 
processes, govf/med both by momentum exchamge. Relevant phenomena are 
expressed by ferms /21-5/ and 721-6/. 

Under evaporative conditions, water flows from the root zone through 
the root epidermis, the' plant hydraulic system and the stomata to 
the atmosphere. In the root zone, water flow meets a resistance depend- 
ing on soil water pressvure. The resistances in root epidermis and plant 
hydraulic system are taken to be constant in this study. The stoma tal 
resistance depends on the opening of the stomata which controls the 
release of water to the atmosphere. 

While water flow is mainly governed by plant physiological factors , 
heat flow is more passive, depending on the plant's ability to evaporate 
'Though the heat flow resistances in soil, canopy and atmosphere are 
variable, they do not influence the heat flow considerably. 

Fig. 5 shows the model's physical diagr 2 un. 

- Parameters and rg 

The model uses a Businger-Dyer approach to estimate the turbalent diffu- 
sion resistance r^ which is given as a function of wind velocity, the 
stability of the atmospheric boundary layer just above the surface, and 
the nature of the surface. Different terms are applied 'under stable, 
neutral and unstable atmospheric conditions. 

Crop resistance can simply take into account stomatal resistance Cg. 

In fact, euticular resistance is at least one order of magnitude higher 
than tg and flow through the cuticula may be neglected. 

Stomatal resistance r^ is a result of closing of the stomata. It is 
caused by decreasing water turgor pressure in the guard-cells surround- 
ing the stomata, which is mainly influenced by leaf water pressure'*/^ 
and-incoming shortwave radiation R . Thus closure of stomata occurs 
when R is low, i.e. during night. 

3 

The stomatal resistance Cg is fo\ind from a set of empirical relations 
with crop height, the short wave radiation and the leaf-water pressure. 


















- The numerical algorithm 


Fig. 6 shows the flow diagram of the algorithm. The soil heat flux is 
estimated using an explicit finite difference method with the spacings 
between each node being equal. The iteration of the surface temperature 
is also simple, the surface temperature being altered by successive 
small amounts until an energy balance has been achieved. This approach 
is satisfactory where the components of the energy balance are not 
changing rapidly between successive time steps, so that the initial 
estimate of the surface temperature (i.e. the temperature estimated 
from the previous time step) is close to the final estimate. 

As noted above, the TERGKA model directly estimates the daily evapora- 
tion. Changes in the soil moisture may also be estimated from the evapo- 
transpiration flux, E, estimated by the term /21-6/. 


- Obtaining soil moisture 

Given the evaporation, soil moisture is obtained by a set of relation - 
ships as follows. For grassland, E may be expressed as 


' E 


i ~ 

plant soil 


725-1/ 


where tp is the leaf water pressure, the soil water pressure, r 

X S Ox3>X11m 

the plant resistance for water transport, and r the soil hydrau-* 

201 1 

lie resistance. Given estimates of , and the resistances r and 

r , . may be estimated. plant 

soil s 


These may be found using empirical relationships with other variables 
and then may be related to tiie volumetric water content, ©, by 



r 


where S is the saturation, defined as 6/6g (6g being the moisture 

content at saturation) . is t.he rest saturation, m is a pore size 
distribution factor, and is a re-scaled soil water pressure; equal 
to where 1/;^ is the air entry value. 
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FIG. 6: FLOW DIAGRAM OF THE TERGRA MODEL 












6.2. 1.3.3. The TELL-US model 


Introduction 

In TELL-US an algorithm has been developed to convert remotely sensed 
day and night temperatures of bare and scarcely vegetated areas into 
parameters relevant to daily evapotranspiration and soil moisture. 

Such an interpretation routine calls for a process which is the reverse 
of temperature simulation. For this purpose, and mainly in the U.S., 
methods have been worked out to calculate the thermal inertia of the 
soil from the difference between the day and the night temperature to 
be obtained by means of HCMM /4, 5 , 6/. This approach has been 
discussed /!/ which is limited. to cases of known or negligible 
small evaporation. Moreover, in this approach only one data point, i.e. 
the D/N temperature difference is used, while two data, the D and N 
temperatures are available. Therefore, in principle, two unknowns could 
be determined. 

This ha,s been achieved by means of the TELL-US algorithm. It uses two 
separate estimates of surface tempera ttire , ideally at the times of 
maximum and minimum thermal 5 tmission (about 14.000 and 02.00 local 
time) and so allows two unknowns to be estimated directly, without the 
need to employ empirical.’ relationships such as those used in TEEGRA to 
derive one of these unknowns. In this model, thermal inertia and evapo- 
ration are estimated; an experimentally-derived relationship may be 
found between thermal inertia and soil moistiire. 

While the general structure of the model is similar to that used in 
TERGEA there are some differences. The soil heat flux is again found 
using an explicit finite difference method, but in TELL-US the spacing 
between nodes is allowed to increase downwards. The iteration at each 
time step is also rather faster than for TEEGRA and it converges rapidly. 

The TELL-US system does contain only two unknowns: the surface relative 
humidity, h, and the thermal intertia p = . 

The algorithm 

The algorithm consists of a main program and a set of sub-programs simu- 
lating the daily course of the surface temperature. The main program com- 
pares_ the simulated D and M temperatures with the remotely sensed temper- 
atures. It produces suitable look-up tables. 

A detailed description of the model is given in / 3 /. 

In Table 10 the essential input data are reviewed. They have been re- 
grouped into "field" and "regional" data. Thefisld, data must be spe- 
cified for every ground plot to be interpreted. When it is deemed valid, 
the regional data may be obtained from the nearest weather station. 


Spatial resolution 
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Table 10 input data for the TELL-US algorithm. 
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Slope direction "and slope dip have to be estimated from topographical 
maps. The surface temperatures and the &lbedo may be obtained by means 
of the HCMM satellite or an airborne scanner. 

Test of the algorithm 

In April 1978, the TSLL-US algorithm has been tested on data from a 
field experiment of the USVICL at Phoenix, Arizona /8 /. Data measured 
2.6 emd 14 days after irrigation have been used. The soil temperature 
measured at 1 mm depth was taken to be a good approximation of the 
surface temperature. The soil moisture profiles were measured gravi- 
metrically. The total daily evaporation was measured by a weighing 
lysimeter. 

Comparisons with the values determined by means of the TELX.-US algorithm 
aure reported in /9 /. The agreement between measured and calculated soil 
moisture is good: it goes from the full coincidence (2 days after irri- 
gation) to an average moisture content of the top of 8 cm of the soil 
(14 days after irrigation). 

Comparison for total daily evaporation is not yet satisfactory as, except 
for the initial drying stage, the algorithm tends to underestimate eva- 
poration. The reason for this tendency lies in the assumption of a 
constant surface relative humidity which causes excessive condensation 
to be simulated for the night. 

6 . 2 . 1 . 3 . 4 Regional evapotranspiration 

The knowledge of actual evapotranspiration (ET) on a regional scale (e.g. 
10 to 10,000 km^) would be of considerable interest in climatology, 
hydrology and agriculture. In particular, its affect on the local 
climate, its role as a major component in the regional water balance 
and its effect on crop yields could be better evaluated. Therefore, 
three approaches specific for the use of HCMM data in estimating regional 
ET have been considered (Seguin, 10/* 

1. Correlating day/night differences in surface temperature with ET 

2. Differential thermal approach. This approach is based on the re- 
lationship between the difference in vapour flux between two sufaces 
and the difference in their surface temperatures, provided the sur- 
faces are similar except for their water regime. Soil heat flux and 
albedo have to be assumed equal for the two surfaces. 

In this case, using equation /2 1.9 /the evapotranspiration 

ET, = ET^ + (4QrT^ + CD h) (T - T ) /29.-1/ 

1 2 a ■ si s2 

from a surface 1, ET^, cam be calculated with the aid of a known 
evapotranspiration from a reference surface, ET 2 , and the difference 
in their surface temperatures, Tg^ - Tg 2 (Itier and Perrier, 11). 
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3. Combined aerodynamic-energy balance approach. This method is 
essentially the same one as used in the TERGRA model. It consists 
in calculating the sensible heat flux by the aerodynamic equation 
relating flux to the difference in air and crop surface temper- 
ature. The latent heat flux* the evapotranspiration ET, is then 
derived from the energy balance. 


6.2.2 Field operations 
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Field operations caurried out in the pre-launch period had the purpose 
of testing the algorithms to be used for the interpretation of HCMM 
data. At the same time, the thermal infrared data produced served to 
prepare the necessary data processing routines and acquainted the 
co-investigators witli the handling and interpretation of this type 
of data. 

Field operations were of two kinds: 

- Flight experiments involving both airborne and ground based measure- 
ments of a duration of days, 

- Continuous field measurements lasting a number of weeks or even 
months. 


6.2.2. 1 Flight experiments 

Three flight experiments took place in the pre-launch phase, in the 
Netherlands, the United Kingdom and in France. While the Dutch flights 
were organized by the National Institutes, the U.K. and French flights 
were Joint Flight Experiments (JFE) with the participation of both 
National Institutes and the Joint Research Centre. 

6. 2. 2. 1.1 The Netherlands (ICW, KNMI) 

Flights with an airborne MSS having 10 channels between 0.38 and l.l^um 
aiid one thermal infrared channel from 8-14 um were performed over two 
areas. The first area, located in the Lopixerwaard (51°38' N, 4*^55' E) 
consisted mainly of grassland on clay soil. The second area was situated 
in the Ijssel Lake Polder (50°13' N, 7°3 ’ E) and was occupied by various 
agricultural crops, especially wheat, sugar beets and potatoes. The 
soil was a silty clay. 

Flight altitude was 1200 m resulting in 2.5x3 m resolution elements on 
the ground. The flights were carried out at the maximum of the daily 
temperature amplitude between 12 and 14 hrs solar time. The dates 


--’f- »• -«»vg 3r , T .,» r »a- f r 
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were May IS, June 14 and August 31, 1977 for Lopikerwaard and May 17 
and June 14, 1977 for the Ijssel Lake Polder. 

Purpose of the flights was testing of the TERGSA model for the evaluation 
of ET and soil moisture from crop surface temperatures. 

Details of the operation, as well as a description of the instrumentation 
on the ground is given in TELLUS Newsletter 2 (Annex XI) . 

Results and conclusions ; On thermal scanner images the relation scan 
angle/sun angle influenced considerably the temperature measured. 

For grassland, as well as for wheat, calculation of evapotranspiration 
from reinotely sensed surface temperatures only was unreliable. Crop (wheat) 
temperature's simulated with the TERGRA model were higher by a few degrees 
them measured temperatures. 

More research should be carried out on the physical interpretation of 
remotely sensed crop temperatures. In the meantime, a thermal differential 
approach (section 3. 2. 1.3. 4) should eventually be used. 


6.2.2. 1.2 France (INRA/ST/SB, JRC) 


A Joint Flight Experiment (JFE/Framce) was carried out on Sept. 30 , 1977 
on the French test site of the Beauce along a flight axis of 20 km 
centered on a point south of Chartres (48®13'49" N, 1*^36'01 E) . The 
flight axis crossed a mixed agricultural area covered mainly by cptn, 
sugar beets, stuibble and bare soil. 


One experimental plot on wheat stubble was equipped for measurements 
of the energy balance. Surface temperatures were meas’.ired on the 
ground with a mobile PRT-5 infrared thermometer. 


Two flights at 1700 m altitude were carried out along the same axis, 
the night flight at 00.52 hrs amd the consecutive day flight at 13.05 hrs 
local time. On-board instirumentation consisted of a 10-channel Daedalus 
MSS, a Daedalus thermal scanner 3-14 .um, a PRT-3 infrared radiometer 


and a 4-channel SAT thermal scanner. 


/ 


The purpose of the JFE was an investigation of the surface temperature 
distribution in various types of vegetation. Details of the operation 
including a description of the instrumentation are given in Annex XI-. 

Results and conclusions ; Surface temperature variations were found to 
be larger inside one crop than between various fields of the same crop. 
During the day, meteorological data needed for the interpretation of 
remotely sensed surface temperatures should be obtain ’5d at exactly the 
same time as the surface temperature. 

A simple model containing three resistances will be used to calculate 
ET. In this model, surface temperatures from HCMM she ild serve to 
correct the values of stoma tal resiS'tance. 


\ 
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Themal data from this flight was used to simulate HCMM resolution 
over a mixed agricultural area. The result of this simulation is re- 
ported in TELLUS Newsletter 1 (Annex I) . noPOnTTCIBlLITY OF THF 

SSai. page is poor 

6. 2. 2. 1.3 united Kingdom (IHW, UL/dg, UR/DG, JRC) 

A Joint Flight Experiment (JFE/UK) was carried out on Sept. 13 on the 
U.K. test site along two flight axes of 10 km each. The first flight 
axis was over the Grendon-Underwood Experimental Catchment (1 W, 

51 53' N) composed of gently sloping land with fields of up to 15 ha. 

The second flight axis passed over the Newbury area (1®44' W, 51®20' N) 
with fields up to 12 ha. 

Soils on the two sites ranged from brown calcareous soils amd surface 
water gley to clays. Lamd use was cropland (cereals) and grassland with ■ 
some forest and market gardening. 

Grassland auid bare soil were instrumented on the Grendon site, 
measuring points on arable soil, grassland and bare soil were estads- 
lished on the Newbury site. A detailed description of the flight and 
of the instrumentation on the ground is given in Annex XI. 

Two consecutive flights at an altitude of 1000 m were carried out along 
the same flight path: A night flight between 04.25 and 05.15 GOT and a 
day flight between 12.45 and 13.35 GMT. On-board instrumentation con- 
sisted of a Daedalus DS-1250 MSS with a thermal channel from 8-14^um 
a Barnes PET-5 infrared thermometer coupled to a Hasselblad camera 
and a RC-883 camera. 

The purpose of the JFE was acquiring experience in the interpretation 
of thermal imagery of agricultural land surfaces. In particular, data 
should be obtained to test the algorithms (TERGRA amd TELL-US) for the 
evaluation of soil moisture amd evapo-transpiration from remotely 
sensed surface temperatures . 

Results auid conclusions ; The thermal and visible-IRR airborne MSS data 
obtained during the campaign were of excellent quality. They were used 
to develop data processing routines li.ke geometric corrections and the 
superposition of day aind night images. They also served as data for 
speeding up the execution of the TELL-US algorithm and for tasting the 
sensitivity of both the TERGRA and TELL-US model. 

An important part of this work has been carried out at the JRC-Ispra 
Image Processing Laboratory which disposes of the following installations: 

- DEC PDP-11/70 with floating point processor and 512 KB core memory, 

- two disk units of 67 MB each, 

- one dual density 1800/1600 3PI tape unit, 

- card reader, line printer, teletype and video terminals, 

- S.E.P. VIZIR high-precision, high-speed laser image recording unit 
(up to 30,000x30,000 points in 64 grey levels with resolutions of 
13.5 u and 40.5 u) , 


- DEC PDF- ll/OS with 56 KS core memory, 

- four logical disk units of 2.5 each, 

- video terminal and line printer, 

- COMTAL image display unit with 3 image memories of 512x512 bytes 
each, contrast function and pseudo-colour tables, trackball /target, 

3 graphic overlay memories. 

The PDP-11 computers are connected through a 1 Mbaud communication 
link. 

TSLL-US model . Soil moisture and evaporation was mapped by processing 
in the TELL-US model pixel by pixel values of day and night surface 
temperature and albedo together with ground data (e.g. surface rough- 
ness) varying field by field and meteorologic data (e.g. solar radiance, 
windspeed, air te;|iperature) valid for the whole flight strip. To this 
purpose, a tabulated version of the algorithm was prepared which uses 
multiple interpolations within tables to reduce processing time to 
a few msec per pixel. The scanner data were first corrected for scan 
amgle effect and for geometry to permit superposition of night thermal 
data on day visible and thermal data. DetadLls of the methods can be 
found in Annex XV . 

Three types of thematic maps were obtained representing thermal inertia, 
soil moisture and cumulative daily evaporation of bare or 
vegetated areas (bare soil, burnt and unbumt stubble, stubble with 
some green shoots) . Each of the three parameters had five to seven 
classes of magnitude and was presented with a ground resolution of 
2. 5x2. 5 m, equal to that of the scanner data. 

It could be confirmed, that a detailed mapping of thermal inertia, 
soil moisture and cumulative daily evaporation was possible using a 
complex model like TELL-US. Obviously, the data processing scheme 
developed is independent of the intrinsic performance of the model 
and each further improvement of the latter will bring some correspond- 
ing improvement of the thematic maps produced. 

The comparison with measured soil moistiu:e values suffers from 
experimental uncertainties, but it seems to indicate a tendency to 
underestimate soil moisture content. Direct measurements of evaporation 
were not available at Grendon, but independent calculations of ET agree 
within 10% with the values obtained by TELL-US. In inhomogeneous areas, 
like burnt stubble, where temperature, albedo and surface roughness vary 
considerably in space, these variations tend partially to mask the 
smoother evolution of moisture and evaporation. 

An analysis of the sensitivity of the TELL-US model to various input 
parameters and to the value of the soil-heat capacity used in the model 
was presented in TELLUS NEWSLETTER 8 (Annex VIII) . calm air, as 
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often found at night, the variations in air temperature due to topo- 
graphy and ground-cover may have an important influence on the thermal 
inertia calculated with the TSLL-US model. 

TERGRA model . A sensitivity amalysis of the model carried out for the 
gzasslauid site at Grendon showed the model not to be sensitive to 
soil parameters at the soil moisture range found on the site 
(Annex III) . 

Cumulative daily evapotreuispiration agreed to within 10% wi^th ET cal- 
culated from Bowen ratio meastorements on the Newbury site. The same 
agreement was found between calculated and measured soil moisture 
on grassland. 

Spectral reflectance . The spectral reflectance in the four Landsat 
bands was measured on the ground on four surfaces during the flight. 

The surfaces were bare soil, stubble, grassland auid burnt stubble. 

These reflectances were compared with the airborne MSS measurements. 

At lower- wavelengths good agreement was found between airborne and 
ground measured reflectance factors. At higher wavelengths, especially 
on grasslamd, there was an increasing discrepancy between the two 
measurements. .A systematic error in the calibration of the airborne 
multispectral scanner seems to be largely responsible for this discre- 
pancy. Details are given in Annex IX. , 

6. 2. 2. 2 Continuous field measurements 

These measurements were carried out to obtain data sets over long 
time periods permitting the comparison between measured evapotranspi- 
ration and evapotranspiration calculated from remotely sensed surface 
temperatures. The results from two measurement-, campaigns, Italy 
(Perugia) and Frsmce (Avignon) are reported. 

6. 2. 2. 2.1 Italy (ISNP, JRC) 

Prom July 13 to September 1, 1977 measurements were carried out on an 
experimental farm near Citti di Gastello, Perugia (43°27' N, 12°15' E) . 
Two grass plots, about 0.75 ha each, were established. One of the plots 
was kept dry, the second one was irrigated by sprinkling. 

The surface temperature and the reflectance in the four LANDSAT bands 
was measxnred from a S-ra tower on each plot using a PRT-5 infrared 
thermometer and an EXOTECH- lOOA radiometer. 

Net radiation, albedo, wind speed and direction, dry and wet bulb 
temperature, air pressure, soil heat flux and soil temperature were 
also recorded at 30-min intervals on bo-th plots. 


The sunmer of 1977 being exceptionally wet in central Italy, no "dry" 
plot could be established. 

Data from one of the plots covering a six day period in August were 
used in the TEHGPA model to calculate evapotranspiration and surface 
temperature. These calculations were compared with measured 


During the day, there was good agreement between simulated and measured 
surface temperatures, with a tendency for simulated crop temperatures to 
be lower tham the measured ones in the morning and higher in the after- 
noon. 

At night, and for low windspeeds, simulated surface temperatures were 
too low and there was no agreement with measured ones. 

These differences being observed at night when evapotranspiration has 
a very small or zero value, no significant effect on the calculated 
evapotranspiration could be observed. 

From a theoretical point of view, the term describing sensible heat 
flux in the TERORA model could be improved by including free con- 
vection, though this modification would be of small practical consequences. 

6. 2. 2, 2. 2 Prance ( INRA/Cib^SS/SE ) 

Measurements similar to the ones in Italy were carried out in the 
period July 15 to September 20, 1977 at an experimental field near 
Avignon-Montfavet (43®55' M, 4°51' E) . A field of short grass (about 
70x200 m) was divided into two plots. One plot was kept dry, the 
second plot was irrigated by sprinkling. 

Global, reflected and atmospheric long-wave radiation as well as wind 
velocity were measured with one set of instruments for both plots. 

At the centre of each plot the following measurements were carried 
out: nfft radiation, soil heat flux, Bowen ratio and surface temperature. 

The surface temperature was measured with a Heiman KT-24 infrared 
thermometer. 

These measurements were recorded with a 1-min scanning interval on a 
data system integrating hourly mean values. 

Data from the period July IS - August 31 were used to evaluate the three 
methods for calculating regional evapotranspiration described in 
Section 6 . 2 . 1 . 3 . 4 . 

- bay/night surface temperature differences, if available every 8 days 
(as 'lost optimistically by HCMM) , could hardly be used to detect 
water availability under the conditions of the test site. This 
statement would be less firm for daily values. 

- No €".-aluation could be made of the differential thermal approach due 
to f.'.e wet summer which suppressed any significant diffei'ences between 
the '-t and dry plots. 
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- The combined aerodynamic-energy balance approach gave a fairly 
good agreement with ET values, hourly as well as daily, calculated 
' with the Bowen ratio method. 

In 1973, the continuous measurements were transferred to dry and wet 
plots in the Crau compatible with HCMM resolution. 

Details of the 1977 campaign and of the plauined activity for 1978 
are given in Annex. X. 


6.3 Heat Budget 

6.3.1 Natural phenomena 

This part of the investigation concerns mesoscale heat budgets not 
particularly affected by man-made heat release. It consists in parti- 
cular in a study of the effect of different land use patterns (open 
fields, fields with hedges) on a regional heat budget, the localization 
of frost prone areas and the correlation of thermal and other cha- 
racteristics of land surfaces and the occurance of hail storms, 

6.3. 1.1 The influence of tree hedges on the surface energy budget (INHA/ST/SAP) 

The present work is carried out by INRA at the French test site in 
Bretagne where significant changes in the mesoclimate have occiurred due 
to the restructuration of land. In this process the characteristic 
tree hedges bordening fields (bocage) were cut amd embankments were 
levelled. 

Two contiguous test sites of 1100 ha each, one- with hedges and one without, 
were investigated by day and night aerial thermography. Areas with 
hedges seem to warm up and cool more rapidly than open areas. 

A simulation of HCMM resolution showed that at night the temperature 
classes identified in airborne images can be fovind back at HCMM reso- 
lution. This is true to a much lesser degree for daytime images. 

The most significamt time for night thermography in this area is before 
midnight, while the night pass of HCMM seems tp be less favourable. 

According to the simulation, thermal inertia was not suited for the 
characterization of the open and hedged areas. 

Details of this study can be found in Annex IV. The rssu its of the 
simulation will be compared wtih HCMM data when this wi,ll be available. 

6.3. 1.2 Hail storms (UCEA, IGA) 

It is extremely difficult to localize hail storms in advance.' At best, 
a forecast and warning can be given for an area of a few millions of 
hectars, while a big storm covers hardly more than 50,000 ha. As a result 
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the installation of preventive systems against hail are extremely 
costly and complex. 

The present research is being carried out by UCSA in coll^lboration 
with IGA in order to correlate the beginning of hail storms with 
thermal, hydrologic and topographic characteristics of land surfaces. 
It is assumed th^t hail storms are formed by spatial temperature dif- 
ferences on the ground, as well as by evaporation from vegetation and 
soil and by atmospheric turbulences (Annex XVI) . 

A attempt is being made: 

1) to spot the birth zones of thunder and hail storms, paying special 
attention to aspects of agricultural planning? 

2) to obtain a localized forecasting of hailstorms. 

Thermal maps obtained by satellite are the only source of information 
available on a mesoscale basis (some 100 km^) to localize the 
"heat island", which are among the main factors originating storms. 
Variations of thermal inertia inducing thermal differences many hours 
before storms are also very important and may bo singled out. 

The HCMM satellite presents favouraJble characteristics for this in- 
vestigation, on condition that suitable time series of thermal soanes 
would be available to co-investigators-s 

The study area chosen by UCSA is within TS No. 3 . The research has 
started but very few suitable scenes have been received so far from 
GSFC. 

Use of TISOS-N and NOAA imagery is also foreseen in this context. 

6. 3. 1.3 Localization of frost prone areas (UF/IG) 

On the German test site near Freiburg work on this subject was carried 
out by UF/IG using airborne night thermal imagery. The investigation 
was carried out on a subregional scale (100 km^) , a local scale (1 km^) 
and on the scale of a single terrace (0.01 km^) . The purpose of the 
investigation was to delimit vineyard sites in frost prone areas. 

On a subregional scale, thermal scanner data from an altitude of 
3600 m permitted to distinguish the vineyard slopes, warmer by 2-3® 
than the cold plain. 

On h local scale, the airborne thermal image.ry was digitally enlarged. 
By mobile measurements of air temperature on 10 clear nights the 
limits of cold surface temperatures and air temperatures in the valley 
bottom could be shown to coincide rather well. 

Flight altitude had to be lowered to 1000-2000 m above ground level 
to obtain results significant on the scale of one vineyard. With, 
thermal data from 750 m hight it was possible to show that large 


terraces, built for the use of mechanized equipment, formed cold 
traps and were unsuited from a micro-climatologic point of view. 

Details of this investigation are given in Annex VI, 


6.3.2 Anthropogenic heat release (UT/IG) 


The long delay in the distribution of HCMM data and the high percentage 
of cloud cover of the images of tdie German test site received, have 
seriously delayed work on this sifsject. 


In the preparatory phase, work was carried out by UF/IG in Germany on 
wind-induced micro-climatic differences detected from thermal scanner 
data. The study treated a mixed urban- rural area. Details are given 
in Annex. V. 




6.4 Working Groups 


Three working groups (WG) were formed in the project to deal with its 
various scientific aspects. Each, of the six groups of Co- Investigators 
is represented by one member in each working group. Additional scientists 
from inside and outside the project are invited to attend the meetings 
of the WG which are held once or twice yearly. 

The working groups cover the following subjects: 


WG I : Planning of flight experiments .This group was responsible for 
elaborating the general outlines of the Joint Plight Experiments. It 
ceased its activity with the execution of the first JFE. 

WG II : Energy balance, soil moisture and evapo- transpiration. Modelling 
and measurements. The group treats the theoretical basis of the project 
and proposes the direction for further development of models and 
instrumentation. It discusses the results of measurement campaigns and is 
active in the planning of such campaigns. 

WG III : Data processing and atmospheric corrections . All questions 
of data handling are treated by this group, such as distribution of 
computer programs, standardization of formats, streamlining of computer models. 
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6.5 Conclusions 

Due to the long delay in the distribution of HCMM data and the rather 
haphazard way these data are- being delivered, the, present report refers 
mainly to work carried out in th'^ pre- launch phase. 

The only interpretation carried out on HCMM images was a visual one. 

Both VIS, pXR and MIR were found to be of good cruality and well suited 
for the identification of geomorphologic features. Only relative gray 
tones were employed in this analysis and no use of absolute temperatures 
was attempted. 

The pre- launch .aircraft program and the measurement campaigns produced 
a number of useful results, especially in the form of algorithms leading 
from thermal images to maps of soil moisture and evapotranspiration . 
Data processing routines such as geometric corrections and superposition 
of night and day aircraft data have been worked out. 

In the next phase of the project the methods developed should be applied 
to the satellite data. Absolute temperature calibration will be of great 
■'importance at this stage as well as the availability of test areas 
compatible with the ground resolution of the satellite. 
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8. ABBREVIATIONS 


AEM Application Explorer Mission 

AL After. Launch (days) . 

ATI Apparent '^ermal Inertia 

CC Cloud Coverage 

CCT Computer Compatiisle Tape 

CETIS Centre. Europ. de Traitjnent de 1 ' Inf ormation Scientif. 

CMS Lannion Centre d' Etude de Meteorologie Spatial de Lannion 
COI Co-Investigator 

DIR Day Infrared 

DV Day Visible 

D Day 

EC European Communities 

ESA European Space Agency 

ET Evapotranspiration 

GMT Greenwich Meridian Time 

GSFC Goddard Space Flight Center. 

GDTA Groupement pour le Developpement de la Tsledet. Aerospatiale 

HCMM Heat Capacity Mapping Mission 

HET HCMM Ehcperimental Team 

HOM Hotine .Oblique Mercator Projection 

HD High Density 

IPP Image Processing, Facility 

HDT High Density Tape 

ID Identification 

IR Infrared 

IFOV Instantaneous Field of View 

IGN Institut Geographique National 

JFE Joint Plight Experiment 

JRC Joint Research Centre 

LAT Latitude . 

LONG Longitude 

MSS Multispeotral Scanner 

NIR Night Infrared 

NT Negative Trsmsparency 

N Night 

NER “ Noise Equivalent Radiance 

ND Night-Day 

PC Priority Coverage 

PI Principal Investigator 

PPL Priority Processing List 

QL Quick Look 

RO Retrospective Order 

SC Spacecraft 
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TSC 

USWCL 

VIS 


Standard Order 
Temperature Difference 
Thermal Inertia 
Track 
Test Site 

Test Site Coordinator 

U.S. Water Conservation Laboratory, Phoenise, AS. 

Visible 

Working Group 




9. SYMBOLS 


0 

«D 

«s 

G 

g 

H 

h 

LE 

m 

P 


S 

Sr 

T 

Ta 

Tc 

t 


Soil heat capacity 

Air heat capacity 

Water vapour pressure in air 

Saturated water vapour pressure 

Soil heat flTix 

Gravitational acceleration 

Sensible heat flux 

Surface relative humidity 

Latent heat flux 

Pore size distribution factor . 

Thermal inertia 

Net radiation . flux 

Aerodynamic resistance 

Stomatal resistance 

Saturation 

Rest saturation 

Temperati^e 

Air temperature 

Crop temperature 

time 

Depth 


y 

s 

5 

d s 

A 

Q 

a 

Vi 

% 

Vt 


Reflection coefficient,, long wave 
Reflection coefficient, short wave 
Psychrometric constant 
Emission , coefficient 
Soil water content 
Saturated soil water content 
Heat conductivity 
Air densi^ 

Ste fan — Boltzmann constant 
Leaf water pressTire 
Soil water pressure 
Re-scaled soil water pressure 
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